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Abstract 

We address the question of whether the upcoming generation of dark matter search experiments 
and colliders will be able to discover if the dark matter in the Universe has two components of 
weakly interacting massive particles (WIMPs). We outline a model- independent approach, and 
we study the specific cases of (1) direct detection with low-background 1 ton noble- gas detectors 
and (2) a 0.5 TeV center of mass energy electron-positron linear collider. We also analyze the 
case of indirect detection via two gamma-ray lines, which would provide a verification of such a 
discovery, although multiple gamma-ray lines can in principle originate from the annihilation of 
a single dark matter particle. For each search "channel", we outline a few assumptions to relate 
the very small set of parameters we consider (defining the masses of the two WIMPs and their 
relative abundance in the overall dark matter density) with the relevant detection rates. We then 
draw general conclusions on which corners of a generic dual-component dark matter scenario can 
be explored with current and next generation experiments. We find that in all channels the ideal 
setup is one where the relative mass splitting between the two WIMP species is of order 1, and 
where the two dark matter components contribute in a ratio close to 1:1 to the overall dark matter 
content of the Universe. Interestingly, in the case of direct detection, future experiments might 
detect multiple states even if only ~10% of the energy-density of dark matter in the Universe is in 
the subdominant species. 
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I. INTRODUCTION 



The microscopic nature of the dominant component of the matter density of the Universe, 
the non-baryonic dark matter, is presently unknown. While several theoretically compelling 
hypotheses have been put forward to explain the fundamental physics of particle dark matter, 
the paradigm of weakly interacting massive particles (WIMPs), has increasingly attracted 
more attention and has been closely scrutinized both from the theoretical and from the 
experimental standpoint (for comprehensive reviews of WIMP dark matter candidates and 
their detection, see, e.g., [U EJ [3j Hj; for a pedagogical review of dark matter astrophysics 
see, e.g., [5]). 

Remarkably, the standard calculation of the thermal freeze-out of species in the early 
Universe [U] predicts that WIMPs generically have a relic density today which is very close 
to the measured density of dark matter ^dm^ 2 — 0.11 [7j. While elusive, WIMPs are in 
principle detectable since in standard scenarios they can interact non-gravitationally with 
ordinary matter via weak interactions. Two approaches have been pursued to detect non- 
graviatational signals from WIMPs in our Galaxy: direct detection and indirect detection. 
Direct detection hinges upon directly observing the scattering of WIMPs off the nuclei of 
ordinary matter in very low background terrestrial environments. Indirect detection, instead, 
aims to detect and identify the annihilation debris of WIMPs which annihilate in the Galaxy, 
including antimatter (positrons, antiprotons, antideuterons), high-energy neutrinos, and 
gamma rays. Last but not least, WIMPs are predicted to exist in several potential extensions 
to the standard model of particle physics including, for example, supersymmetry, models 
with universal extra dimensions or little Higgs models. The new physics expected in these 
models can manifest itself as distinctive signatures at high-energy colliders and WIMPs, 
in many cases, are the final result of (potentially long) decay chains when new particle 
states are pair-produced in high-energy collisions. A typical signature at colliders includes a 
large missing transverse momentum associated with the WIMP escaping the detector, and 
a characteristic series of energetic jets and leptons streaming off of the WIMP-generating 
decay chain. 

When exploring the realm of unknown physics it is certainly a reasonable approach to 
invoke Occam's razor and to assume as a working hypothesis the simplest scenario consistent 
with observations. Given the persuasive theoretical and observational indications discussed 
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above, if we restrict our attention to a WIMP dark matter scenario, it makes sense to 
postulate that a single stable WIMP species accounts for all of the dark matter. Yet, while 
a perfectly acceptable and reasonable assumption, this might be an oversimplification of 
the truth. The dark-matter sector might not be as simple as to be made up by one single 
component, and assuming so may obscure new physics in the dark sector. We believe it 
is legitimate to entertain the hypothesis that dark-sector physics reflects the variety and 
complexity of elementary constituents that populate the ordinary-matter sector, and it is 
this hypothesis we scrutinize in this work. 

Scenarios with more than one dark matter particle constituent have been envisioned and 



investigated in the past (see our review and discussion in sec. VI). In this paper we focus on 



the possibility that the dark matter has two WIMP constituents. We have two goals: 

• we want to establish whether current experiments that focus on WIMP detection have 
the capability to actually discriminate more than one WIMP species; 

• we want to answer the phenomenological question: if signatures of more than one 
WIMP species can be detected, which portions of the WIMP parameter space are 
favored? 

In the analysis we keep in mind several theoretical principles that could guide model building 
of a multiple WIMP scenario, where the WIMPs have similar physical properties. 

On more quantitative grounds, the scope of the present study is twofold. Firstly, we aim to 
identify which multi-particle dark matter models might be discovered with present and next 
generation dark matter detection and collider experiments; secondly, we aim to compare and 
cross-correlate — using appropriate assumptions to relate the relevant physical parameters 
- the potential of direct, indirect and collider searches for multi-component particle dark 
matter. We must stress that, whereas direct detection and collider searches could lead to 
the discovery of a dual- component scenario, indirect detection could not provide a discovery 
in its own right. This has to do with a fundamental problem related to the detection of 
multiple monochromatic lines: given two lines, one cannot distinguish a second line arising 
from a second WIMP from an additional line due to other annihilation channels of the first 
WIMP (annihilations into Z7 or i/7 for example). But if the discovery is made through 
the other search channels (direct detection and collider), then indirect detection via two 
gamma-ray lines could provide a valuable verification. Keeping in mind this caveat, in this 



paper we also analyze the indirect detection case in detail, in order to compare to the direct 
detection and collider results. 

The manuscript is organized as follows: In Section [TT] we describe the physical parameters 
of our model and survey our main results. SectionlTTllfocuses on direct dark matter detection, 



Section |IV| on collider searches, Section |V] on indirect detection. In Section |VT] we comment 

we 



on related work and potential multi-particle dark matter models, and in Section VII 
conclude with a summary. 



II. MODEL SETUP AND SUMMARY OF RESULTS 

We consider a scenario where two particle species xi an d Xi are stable on cosmological 
scales, and contribute to the total dark matter energy density Odm with energy densities, 
respectively, Q X1 and X2 . We assume that the whole of the non-baryonic dark matter is 
made up of Xi an d X2 particles, i. e. Odm = Q X i + ^X2 j an d that both particles are weakly 
interacting massive particles. We define the ratio of relic abundances of the two species as 

(1) 

"X2 

which allows us write 

w I 
fl xi = — — flr>M and X2 = - — ■ — DM - (2) 
1 + w 1 + w 

As long as (cxi,2Xi,2 u )tot is dominated by its s-wave amplitude we have the well-known 
relation [6]: 

(^1,2x1,2 «)tot(-T = 0) 

which is independent of mass up to logarithmic corrections. The parameter w is thus also 
typically equal to to the ratio of total annihilation cross sections r as 

(^X2X2^)tot(T = 0) 



r 



w. (4) 



Kixi^)tot(T = 0) 

For most of the following discussion we take particle 1 to be lighter than particle 2, i.e. 
ui2 > mi and we define a parameter A corresponding to the relative mass splitting 

A ee m2 ~ mi . (5) 
mi 

In some cases we will trade A for m2, relax the assumption m 2 > mi, and plot m 2 vs mi. 



Given this setup, we wish to address the question of the detectability of at least one of 
the dark matter species, and the question of our ability to discriminate between a single 
and a dual-component dark matter scenario, in the case that both species are detectable. 
In order to proceed, we shall need to outline sets of minimal, further assumptions, pertinent 
to various dark matter detection techniques, as well as to collider searches. In particular, in 
this study we will address: 

(I) the direct detection of dark matter, proceeding via the detection of recoil energy re- 
sulting from the elastic scattering of WIMPs off of nucleons (for a review see e.g. 

0); 

(II) collider searches for new physics, featuring the end-production of two stable, weakly 
interacting massive particles that escape detection carrying away significant missing 
energy and momentum, with a high center of mass energy e + e~ linear collider. For 
simplicity, we do not consider here the otherwise extremely interesting phenomenology 
of scenarios with two stable weakly interacting states that can escape detection at the 
LHC, since discussing in any detail that scenario would require a number of highly 
model-dependent assumptions (including the detailed mass spectrum of particles that 
participate in the cascade decays of the produced primary particle pairs, branching 
fractions, decay modes, etc.); 

(Ill) the indirect detection of dark matter, specifically via the (likely) cleanest possible signa- 
ture of WIMP pair annihilation, namely the detection of a monochromatic gamma-ray 
line from WIMP pair annihilation into two photons (with energies equal to the WIMP 
mass) with the Fermi Large Area Telescope (LAT) [9]. 

We define below the additional parameters needed, within the minimal and model- 
independent particle setup we utilize here, to quantitatively discuss the three experimental 
probes of dark-matter physics listed above. 

A. Direct Detection 

The rate at which WIMPs elastically transfer energy to nuclei depends on astrophysical 
parameters, as the local density and velocity distribution of WIMPs, and on particle physics 
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parameters, as the WIMP mass and its scattering cross section off protons and neutrons. 
We assume here, for simplicity, that the scattering cross section of either dark matter par- 
ticle species is identical for both neutrons and protons 2n ), and consider the following 
parameter space: 

Direct detection Parameter Space : mi, A or m%, w, cr Xin , ^xm- (6) 

Given that mi and A specify the WIMPs masses, and w the relative abundance (which we 
assume to be the same locally and on cosmological scales), the parameter space above fully 
specifies the WIMP properties relevant to calculate the direct detection rates. 

In particular cases, it can be useful to postulate a relationship between the WIMPs 
pair annihilation cross section and the scattering cross section of WIMPs off nuclei, and 
between the scattering cross sections of the two WIMP species under consideration here. 
As a guideline, one can assume that the particles exchanged in WIMP-nucleon scattering 
are the same ones that contribute (in a crossed-channel) to the WIMP pair annihilation in 
the Early Universe. With the case of supersymmetry in mind, if we assume this exchanged 
particle is a "squark" % = 1, 2, and that the squark mass m^ is much larger than the DM 



mass, we can write: 



and 
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a xiXi ~ — r- ( 8 ) 



One therefore can estimate 

r xn = ~ f X2X2 ~ -^w. 9 
°xm m 2 X2 (a XlX1 v) m? X2 

When appropriate and explicitly specified we will assume the above relation to hold 
and use it to relate one scattering cross section to the other according to Eq. ([9]), e.g.: 

a X2n ~ wa Xin {m x Jm X2 f. 



B. Collider Searches 

Outlining the phenomenology of models with two stable massive weakly interacting particles 
that end the decay chain of new particle species produced in high-energy collisions is a 
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formidable and complex task. Even simply counting event rates that involve one or the other 
particle involves specifying a number of features of the particle model at hand, including 
the production cross section for the heavy particles that subsequently decay into the stable 
WIMPs and the topology and branching ratios for the decay chains that lead to the escaping 
WIMP(s). In addition, in the present analysis we aim at understanding whether it is possible 
to not only detect the existence of WIMPs, but to discover more than one stable WIMP 
species and to discriminate this scenario from one where there is only a single stable WIMP. 
For a model-independent analysis, in the case of colliders, an unambiguous experimental 
handle on the particle mass is likely needed. Although several studies argue that, for specific 
models or event topologies the Large Hadron Collider (LHC) will provide such handle in 
single WIMP models (see e.g. [TU] for extensive discussions in the case of supersymmetric 
models, and Ref. [Uj for specific event topologies), it is obvious that a hadron collider is not 
the ideal environment for precision studies of the absolute mass scale of new particles. The 
main reason is simply that the center-of-mass energy of the underlying event in a hadron 
collider is not known. We leave the question of model-specific scenarios with multiple stable 
WIMPs at the LHC for future dedicated studies. 

We instead articulate in the present study an analysis of a high-energy electron-positron 
collider that might follow the LHC at some point in the future. With such a machine it will in 
general be possible to tune the center of mass energy of collisions, and to have an extremely 
accurate account of the mass scale of new particles. In this case, the ability to statistically 
detect a signal and to discriminate between a single and multiple stable WIMP scenarios 
(where stable naturally alludes to collider time-scales) is a well defined question that depends 
only on the number of events produced in high-energy collisions and assumptions about the 
experimental features of the collider detectors and relevant backgrounds. We give details on 



our assumptions about all these parameters in Sec. |IV[ where we consider, as a prototypical 
future electron-positron collider, the so-called International Linear Collider (ILC). 



While we give full details of our assumptions in Sec. |IV[ we wish to first briefly outline 
here the scenario and parameter space we shall consider. As it is unlikely that stable WIMPs 
will be the most easily produced particles at the collider, we envision a simple scenario where 
a pair of heavy charged particles C of mass are produced that subsequently decay into 
either xi or X2 particles. In the limit that processes involving the C particle dominate 
the total pair annihilation cross section (for instance via t— or u— channel annihilation 
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into standard model fermions, as is the case of neutralino pair -annihilation into fermion- 
antifermion final states via sfermion exchange) we have the relation 

r a2 = 4-£^4^ (10) 

between the ratio of the squared structure constants r a 2 and the ratio of relic abundances w. 
The relevant Feynman diagrams are shown in Figure u\ Here 0:1,2 = 91,2 is the coupling 





Xl,2 



(a) 



(b) 



(c) 



FIG. 1: Feynman diagrams. 1(a) Pair production of at collider, 1(b) Dark matter annihilation, 



1(c) Decay of C ± . 



constant for the vertex involving C ± , X12, c*, where denotes a charged standard model 



particle. Using Eq. (10), we can then relate the branching ratios for C decay to xi,2 to w. 



Accounting for the two WIMP masses (determined by mi and the relative mass splitting 
A), the simplest parameter space we can envision for this scenario is: 



Linear Collider Parameter Space : m 1; A or m 2 , r a 2 ~ w, mg. 



(11) 



C. Indirect Detection: Monoenergetic 7 Rays 



In the era of precision gamma-ray astronomy, a golden indirect detection channel is the 
search for monochromatic 7 rays resulting from the pair annihilation of two dark matter 
particles into two photons (for early studies see, e.g., Ref. [12], see Ref. [13] for the complete 
calculation in the context of the minimal supersymmetric extension of the Standard Model 
(MSSM), and Ref. [H] for the extension of that calculation to the next-to-MSSM). For non- 
relativistic (cold) dark matter particles, which pair-annihilate at very low relative velocities, 
the energy of the two resulting 7 rays equals the mass of the dark matter particle. 
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While in specific WIMP models the monochromatic gamma-ray line is a suppressed an- 
nihilation mode, it is by all means a very striking signature that could not be attributed to 
astrophysical processes. We might have considered the gamma-ray continuum emission as 
well, which typically has large fluxes. It is clear, however, that discriminating a dark matter 
signature from more than one WIMP from astrophysical backgrounds (be it from data taken 
from the galactic center or from the galactic halo diffuse emission) will be a hopeless task for 
an instrument like Fermi. Since we want to be here agnostic about other pair-annihilation 
channels, that would be needed to compute the continuum gamma-ray emission, and in view 
of the difficulties of background discrimination, we focus only on the monochromatic line. 

We can still not claim any discovery of a dual- component dark matter via detection of 
multiple lines, because of the complication discussed in the introduction, but this remains 
the cleanest signal that we could use as supporting evidence if such a scenario is found via 
other searches. 

We could have also considered another smoking gun signature of dark matter annihilation: 
a flux of high-energy neutrinos from the Sun or from the center of the Earth. Again, though, 
this would imply the need for (a) detailed assumptions on the WIMP annihilation modes, and 
thus would introduce model-dependence, and (b) it is unclear if the spectral information one 
gets from current or future neutrino telescopes would ever allow for a robust discrimination 
of scenarios with more than one WIMP. Though interesting, we thus choose not to entertain 
this indirect detection channel in the present analysis. 

Given that dark matter is dark (i.e. it has suppressed electro- magnetic interactions), the 
pair annihilation cross section of WIMPs into two monochromatic photons is, in generic 
WIMP models PQ [21 El S], a loop-suppressed process. In specific constructions, however, 
the monochromatic photon pair-production process can be enhanced, as shown for instance 
in the case of inert Higgs dark matter scenarios [15], models with extended gauge sectors 
and higher dimensional operators [16] and string-inspired scenarios [17]. In the present 
model-independent approach, we indicate the relevant annihilation cross sections with the 
symbols: 

o-(Xi,2 Xi,2 ^ll) = 0- xli2Xli2 ^ 7 • (12) 

For simplicity we assume that the process where X1X2 77 is suppressed or forbidden and 
we neglect it (given this would give an additional line signature, this is a conservative as- 
sumption). To obtain a more tractable parameter space we reduce the number of parameters 
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and make the simplifying assumption that the total rate of pair-annihilation is proportional 
to the rate at which the WIMPs pair annihilate into two photons, 



a Xi, 2X1,2—77 ^ ( a Xi,2Xi,2 V )tot(T — 0). (13) 

This assumption is valid in the context of several particle dark matter models: for instance, 
in models such as those with universal extra dimensions where the new physics sector lies at 
a quasi degenerate mass scale tokk, both (cr xx t>) to t and <r xx ^ 77 are proportional to J7% K [3]. 
Similarly, Ref. (18] showed that in the case of supersymmetric Higgsinos (another paradig- 
matic "minimal dark matter" particle model [IS]), the ratio of the two cross sections is 
constant, modulo logarithmic corrections in the ratio of the weak scale to the particle dark 
matter mass scale. The approximate result of [18J was confirmed by the exact calculation of 
Ref. [20J. In addition, several examples of model-independent constructions where the two- 
photon mode has a roughly constant branching ratio to the total annihilation cross section 
are given in Ref. |21j . 

Using Eq. Q, we can then argue that the ratio of the annihilation cross sections to two 
photons is 



r 7J = ax2X2 ^ ~w. (14) 

^XlXl— 77 

This set of assumptions is generically satisfied in most particle dark matter models, 
unless special circumstances occur, e.g. coannihilation processes, or accidental suppressions 
or enhancements of the radiative pair-annihilation channel into two photons. Once the mass 
of the WIMPs is specified, and the pair annihilation rate into two photons is given, the 
rate for the indirect detection channel under consideration here is fully specified from the 
particle physics standpoint. We therefore take the parameter space for this indirect detection 
channel to be: 

Indirect detection Parameter Space : mi, A or 7712, r 77 ~ w, (o" XlXl ^ 77 t>). (15) 

Our ability to make predictions therefore only depends, once these parameters are known, 
upon specifying the dark matter density distribution in a given direction, the astrophysi- 
cal background, and the specifications of the instrumental response and performance. We 
provide details on all these, and present our results on this parameter space in Sec. |Vl 
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FIG. 2: Summary of our results: A highlight of the results we present in this study for direct 
detection (upper panels), collider searches (central panels), indirect dark matter detection (lower 
panels). Color coding: blue stands for no detection, green for detection without discrimination, red 
and yellow for detection with discrimination. 

D. Summary of the Results 



For the ease of the reader we now briefly summarize our findings and outline the main 
highlights of our results in this subsection. Complete details on indirect detection, direct 
detection and collider searches are then given in the following three sections of the paper. 
Figure [2] shows a selection of representative plots from the main body of the paper. We 
also take this opportunity to introduce the color-coding and conventions we shall use in the 
remainder of this study. 
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The top panels of Fig. [2] refer to direct detection (the results refer to a prototypical 
experiment employing 1 ton of Xenon and 5 years of data taking time), the middle panels 
to collider searches, the bottom panels to indirect dark matter searches. In all cases, the 
panels to the left refer to a WIMP setup where A = 1, hence the stable WIMP masses are in 
a 2:1 ratio (m 2 /m 1 = 2), and where they contribute equally to the dark matter abundance 
w = Q X1 /Q X2 = 1; the parameter space refer to the mass scale {m\ = m 2 /2) for the x-axis, 
and to the variable that drives the event rate on the y-axis: respectively, the scattering rate 
off of nucleons for direct detection, the mass of the heavy species to be pair produced at 
a linear collider for collider searches, and the pair annihilation rate into two photons for 
indirect detection. The right panels inspect alternative slices of parameter space for each 
detection technique. 

The color coding we employ here and in the rest of the paper is as follows: 

• Red regions correspond to detection with discrimination: a combination of parameters 
such that we predict that the experiments should both detect a signal from both WIMP 
species, and discriminate the individual WIMPs (for instance, by distinguishing the 
two different WIMP masses in direct detection and collider searches). 

• Green regions correspond to detection without discrimination: in these regions of the 
parameter space, at least one of the WIMPs is conclusively discovered, but it is not 
possible to disentangle the presence of a second WIMP species. This is either be- 
cause the second WIMP doesn't produce a large enough signal, or because the signal 
it produces cannot be experimentally distinguished from that coming from the first 
WIMP. 

• Blue regions indicate portions of the parameter space with no detection, i.e. where 
the experimental setups we consider would not see any statistically significant signal 
from either WIMP species. 

In the case of direct detection, assumptions on the systematic control of background play 
an important role. We therefore also consider a 

• Yellow region, where one would have detection with discrimination (as in the red 
regions) in the optimistic limit of perfect background control, i.e. if one assumes that 
the rate of background events is perfectly known. 
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This clearly extends the red region to wider regions of parameter space, where an imperfect 
background knowledge would impair the experimental ability to discriminate the existence 
of two WIMP species. 

In the top left panel we learn that, compatibly with current direct detection limits in 
the parameter space of mass versus scattering cross section off nucleons, future detectors 
might be able to conclusively detect two WIMPs making up half of the cosmological dark 
matter each, as long as the mass of the lighter WIMP is below 40-50 GeV, and the heavier 
WIMP is lighter than 100 GeV or so. Also, the scattering cross section must not be much 
beyond current limits. The yellow shaded area indicates that while a good control on the 
background will not extend the discovery reach for a multi-WIMP scenario in scattering 
rate, it will indeed enlarge the detectable parameter space in mass, when current limits are 
considered. The green region in the figure points out that to detect two WIMPs might 
require scattering rates which are much larger (by roughly two orders of magnitude) than 
those needed to simply discover one of the dark matter particles. The panel to the right 
highlights, for a fixed value of the cross section, the regions, in the (mi, 1712) parameter space, 
where discovery of a multi-particle dark matter sector is possible (the grey shaded area is 
excluded by current limits). There exist two general regions of parameter space where one 
can expect to discover and disentangle two WIMPs with a ton-sized Xenon detector: one 
corresponds to the lighter species having a mass around 10-20 GeV, and the second one in 
the 30-50 GeV range. The second region features a light mass for m 2 and extends in a wide 
range of masses for mi (which is assumed to have a relatively large scattering cross section, 
10~ 44 cm 2 ). In short, a general requirement is that one of the two WIMPs be light (tens of 
GeV), and the second one sufficiently strongly interacting with matter. 

The central panels explore the capabilities of a linear collider with a center of mass energy 
of 0.5 TeV. The contours reflect pretty obviously the requirements of kinematics: The blue 
region corresponds roughly to where mc < mi, and thus C cannot decay into either WIMP. 
In the green range, we approximately have mc < 2 x mi = m 2 , so that only mi is produced. 
When m c > mi, m 2 essentially the entire parameter space can be explored with a linear 
collider and the two stable WIMPs unambiguously discovered. In the right panel we point 
out that the only caveat to this situation is when either the two masses are very close to 
each other (i.e. small A) or they are close to either mc or the kinematic reach of the collider 
(respectively, the upper left portion of the red contour, and the right edge). Still, a linear 
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collider would be potentially able to discover two stable WIMPs even if the relative mass 
splitting is of the order of a few percent (left edge of the red region). 

For indirect detection, we indicate with dashed lines in the two lower panels the lim- 
its on monochromatic gamma-ray lines from current, pre- Fermi data, namely EGRET (at 
low masses) and atmospheric Cherenkov Telescopes (larger masses) [221 122]. We will give 
details of our assumptions in the relative section. The left panel shows that in this chan- 
nel the prospects for singling out two gamma-ray lines originating from the monochromatic 
Xi, 2X1,2 - * 77 channels are rather promising. The panel to the right highlights the existence 
of the "sweet spot" for Fermi-LAT in the parameter space of the pair annihilation cross 
sections, setting the first mass to 40 and the second to 80 GeV. For cross sections of the 
order of 10~ 25 cm 3 /s, and this combination of WIMP masses, Fermi-LAT can disentangle the 
existence of two monochromatic lines, with cross sections compatible with current gamma- 
ray data. The two green regions correspond to where we expect to detect either the lighter 
WIMP xi (upper green region) or the heavier WIMP \2 (green region to the right). 

III. DIRECT DETECTION 

A large class of current direct dark matter detection experiments measures the number 
iV of elastic collisions between WIMPs and target nuclei in a detector, per unit detector 
mass and per unit time, as a function of the nuclear recoil energy E r . Given a local dark 
matter density po — 0.3 GeV cm -3 (here with some fraction due to the first WIMP and 
another fraction due to the second one) and the velocity distribution f{v x ) of WIMPs near 
the earth, the differential rate per unit detector mass and per unit time can be written as 



where the WIMP-nucleus cross section, <j x n, is related to the WIMP-nucleon cross section, 
cr x „, by a X N = a xn (Am r /M r ) 2 , with M r = ™ x ™^ the WIMP-nucleon reduced mass, m r = 
mx , mjv the WIMP-nucleus reduced mass, m v the WIMP mass, m„ the nucleon mass, tun 
the nucleus mass, and A the atomic weight. F(E r ) is the Woods-Saxon form factor [23] 



[24j 




(16) 



F(E r ) 



3ji(ggi) l (gs) 2 
qRi 



(17) 
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where q = y/2m^E r is the momentum transferred, ji is a spherical Bessel function, Ri = 
VR 2 - 5s 2 , with R ~ 1.2 A 1 / 3 fm, A the mass number, and s ~ 1 fm. For the velocity 
distibution we take 

/K)^x = (18) 

where u° ~ 220 km/s is the velocity of the Sun around the galactic center, and we have 
neglected the motion of the Earth around the Sun. After integrating over the angular part 
in order to find the speed distribution one gets 

4t> 2 

fMdv x = T-^ e ~ {VxK?d ^ (19) 



The integration over velocities in Eq.(16) is limited, for a given recoil energy, to those that are 



kinematically allowed, meaning that there is a minimal velocity given by v m i n (E r ) = y m 2 ^f r ■ 
From Ref. [23], we take the background to be at the fixed value of N bkg = 1.2 x 10~ 4 
events/kg/day. 

To asses whether a WIMP setup is or not detectable and whether we can discriminate 
between a single and a dual-component dark matter scenario, we proceed in the following 



way: we compute the differential scattering event rate (16) for a given experiment and data 
acquisition time, for the two-component DM setup with the parameters of Eq.([6]). We then 
use a x 2 criterion to tell one of the three situations apart: (1) no detection: the fit to 
the background is better than the given confidence level; (2) detection, but no discrimina- 
tion: the fit to the background is worse than the given confidence level; (3) detection and 
discrimination: the single best fit DM model does not give a good enough fit. 

Let us explain in more detail how the \ 2 criterion is used here. N sign is the signal, N bkg 
is the background, and N obs = N sign + N bkg is the total signal measured by the detector. 
We divide the recoil energy range between 4 and 32 keV in n = 7 equidistant energy bins. 
We then compute three different \ 2 '- 

1. Fit using only the background with one parameter. 



If the fit is better than a 95% confidence level we declare no detection, otherwise we 
declare detection. 
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2. Fit using a single DM model with two parameters. 

(fi(a, m) + N bk 9 - N° bs ) 2 



X 2 (o-,m) = 



jyobs 



1=1 

Here, the parameter a is the WIMP-nucleus cross section, and fi(a,m) is the differ- 



ential rate in Eq.(16) integrated over the bin width, with m the mass of the WIMP. 
Note that the background N bkg is held fixed. If this fit is better than a 95% confidence 
level, we conclude that we cannot distinguish between the two species of dark matter, 
otherwise we declare discrimination. 

3. Fit using a single DM model with three parameters. 

i=l » 

The difference between this x 2 an d the previous one is that the background has now 
become a parameter that we are going to fit. Again, we declare discrimination if this 
fit is worse than a 95% confidence level. 

To help visualize how we distiguish among (1) no detection, (2) detection, (3) discrimi- 
nation, we plot in Fig. [3] the number of events measured, for different WIMP-nucleus cross 
sections, when the masses of Xi an d Xi are very close. 

In our study we consider a XENON-type experiment with 1000 kg of xenon and 5 years of 
data acquisition. In the following plots, the colors yellow and red both denote discrimination. 
In the yellow regions we use a chi-square test with two parameters, which corresponds to 
say that one has perfect control over the background; in the red regions we use a chi-square 
test with three parameters, one of which is the background itself. 

Fig. [4] sets the stage showing the parameter space of the lighter WIMP mass, m 1; versus 
the WIMP-nucleon scattering cross section c Xin , for selected values of the second WIMP 
mass (we remind the reader that m 2 = mi x (1 + A)) and of the scattering cross section 
off of nucleons for the second WIMP, set by the relation a X2Tl = w x <j Xlfl (mi/m 2 ) 2 , as in 
Eq. Q. The various pair of values for A and w we employ is specified in each panel, and 
we indicate the current limits with a dashed black line (the region excluded lies above the 
line) . 

We see that for direct detection, a necessary condition for discrimination (whether or not 
the background is well known) is to have a significant mass splitting between the WIMPs. 
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FIG. 3: The number of events per kg per day is plotted as a function of the nuclear recoil energy 
on a logarithmic scale. For a small cross section (dashed lines), both signals for xi arj -d X2 are 
below the background and there is no detection. For a cross section of 10 2 cm 2 (diamonds) we 
can detect a signal, but the two lines, red and black, are too close to each other and we cannot 
distinguish between the two particles. For a large cross section (circles) we are way above the 
background, and now the two lines are far enough from each other to allow for discrimination. 

For A = 0.1 (two top panels), discrimination is 2-3 orders of mangitude above the cur- 
rent exclusion limits. For A = 1, discrimination is possible with a good knowledge of the 
background and for light WIMPs, both in the 10-40 GeV range, with large scattering cross 
sections, around the current limits. Larger mass splittings (A = 10, third line of panels) 
entail detectability with discrimination only if the lighter WIMP is around m,\ ~ 10 GeV 
and the two dark matter candidates contribute comparably to the dark matter content of the 
Universe. Interestingly, since for direct dark matter detection the relevant particle density 
enters the detection rate linearly and not quadratically (as for the annihilation rate relevant 
for indirect detection), large values of w (i.e. a setup where dark matter is dominated by the 
light component) can lead to discrimination of a multi-partite dark matter scenario. For a 
xenon detector, the best mass range for this case lies in the tens of GeV range for the lighter 
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FIG. 4: Direct detection. o~ xin vs mi is plotted at fixed values of A and w. The current limits 
from XENON10 2007 [26] are also shown, with the arrows indicating the excluded region. Color 
coding: blue stands for no detection, green for detection without discrimination, red and yellow for 
detection with discrimination. 



dark matter particle (larger values being favored with smaller mass splittings between the 
two WIMPs, see the bottom panels). 




FIG. 5: Direct detection. o~ xin vs A is plotted at fixed values of m\ and w. The current limits 
from XENON10 2007 [26] are also shown, with the arrows indicating the excluded region. Color 
coding: blue stands for no detection, green for detection without discrimination, red and yellow for 
detection with discrimination. 

In Fig. [5] we explore an alternative slice of parameter space, and investigate the prospects 
for discrimination on the (A, er Xin ) plane, for fixed values of m\ = 40 GeV (upper panels) 
and 100 GeV (lower panels). The limits from current searches are also shown, and lie above 
the region where discrimination is possible only marginally, in the m\ = 40 GeV case, for 
large w = 10, when A is close to 1. In general, we find that the favored range for A for 
discrimination is always very peaked around 0.5-2, meaning that both a much heavier and a 
quasi-degenerate double-WIMP setup will be extremely hard to discover at direct detection: 



19 



-38 



10 " gi m^ mil] imi| mil] mil] imii ] nun] hi 

_.n e m = 40 GeV, A = 0.1 
1Q -4U ^ - 



lo" 42 1 



c _ -44 _ 

** 10 r 



10 



10 



10 



-46 



48 



38 



DIRECT DETECTION 

T = 5 years, M = 1000 kg 

fT iqUliq Hill] llllll] 11111]! I III] II, 



XENON 10 2007 (Nel 136 kg-d) 



=JL__jL_[__i____jJ= =Lt____l 



' excluded 



I I ■ ■ I 



I i i i III" 



m = 40 GeV, A = 1 



nu n] i niq mil] nun] nil] i mn 



| excluded 



ft Sit 

1 



I ■ ■ ' ' I 



mi 1 I 



m„ =40 GeV, A= 10 



^ "I 



detection 



i I "' 



-46 -44 -42 -40 

10 10 10 10 

a (cm ) 



-46 -44 -42 -40 

10 10 10 10 

a (cm ) 

X 2 n 



-46 -44 -42 -40 

10 10 10 10 

a (cm ) 



|J i 1 i ml 1Q" 48 



10 



10 



10 



1 10 



1 10 



■38 



-40 



-42 



-44 



-46 



lO" 40 ki 



1 io- 42 H 



c „ -44 _ 

* 10 r 



10 
10 



-46 



■48 



m n] mil] imi| inn] i mil] ii mil 1 



m = 100 GeV, A = 0.1 



HUH] Hill] Hill] III! 



I i J J limn 



III H] 1 111^ Hill] Hill] Hill] mm mn i mil i inn ; 



m = 100 GeV, A= 1 



111^ llllllj I 



excluded I 



' I mill I i 



nJ I i ii" 



mn] i ni^ mil] i mi] nun] nun 



m =100 GeV, A =10 



llll] llllllj I 



10 
10 

ii 10 



nJ J J I i i 



1 10 
"I 10 



-46 -44 -42 -40 

10 10 10 10 

2 

a (cm ) 



-46 -44 -42 -40 

10 10 10 10 

a (cm ) 



-46 -44 -42 -40 

10 10 10 10 

a (cm ) 



10 



38 



-40 



-42 



-44 



-46 



-48 



FIG. 6: Direct detection. Here we relax the assumption (|9| and we plot <J Xin vs o X2n at fixed 
values of m\ and A. The current limits from XENON10 2007 [26] are also shown, with the arrows 
indicating the excluded region. Color coding: blue stands for no detection, green for detection 
without discrimination, red and yellow for detection with discrimination. 



the best case scenario is for a mass ratio for the two WIMPs close to 2:1. 

In Fig. [6] we relax the assumption of Eq. (|9| and study the ranges of WIMP-nucleon 
scattering cross sections that can be detected, for fixed values of the masses. In the upper 
panels the lighter WIMP weighs 40 GeV and the second one has a relative mass splitting 
of A = 0.1, 1 and 10, left to right. The second row has the same splittings, but for a light 
WIMP of 100 GeV. The panels highlight that the most favorable situation is either if one of 
the two WIMPs has a very large scattering rate compared to the second (the vertical bands 
to the right and at the top of all panels), or if the two scattering cross sections are very 
close to each other and large enough (the wedge in the upper right corners). The size of 
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FIG. 7: Direct detection. Here we replace A with 777-2 m the parameter space ([6]) and plot 777,2 vs tti\ 
at fixed values of o~ Xin and w. The shaded area is the region excluded by the limits from XENON10 
2007 [26 . Color coding: blue stands for no detection, green for detection without discrimination, 
red and yellow for detection with discrimination. 

the yellow regions highlights the importance of background suppression and control for the 
discovery of a multi-partite dark matter setup. 

As evident from Figs. [5] and [6j the viability of a multi-component dark matter scenario 
discovery at direct detection experiments depends very sensitively on the choice of masses, 
and although it also depends on the relative mass splitting and relative contribution to the 
dark matter density, the precise value of the WIMP masses appears to be a crucial factor. 
We explore this fact in Fig. [7j where we plot ni2 vs rri\ (we shade in grey the regions ruled 
out by current experiments [26]). We take a Xin = 10~ 44 cm 2 , and vary w = 0.1, 1, 10 and 
100. 

We find that discrimination is possible if one of the two WIMPs is light, namely with 
a mass below 50 GeV, and the second WIMP has a mass between 20 GeV and 1 TeV. 
These ranges depend on the relative abundance w, but the maximal extension of the regions 
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T = 5 years, M = 1000 kg 




AAA 

FIG. 8: Direct detection. We compare three experiments running for the same exposure time 
and using the same mass of xenon, germanium and argon respectively. Color coding: blue stands 
for no detection, green for detection without discrimination, red and yellow for detection with 
discrimination. 

outside the current experimental bounds is found for large w (i.e. for a Universe dominated, 
in terms of relative matter abundance, by the lighter WIMP). Notice that the discovery 
regions in yellow and red are not symmetric in the plot because of the assumed scaling of 
a X2n with w and the WIMP masses. 

Fig.|8]expl ores the different prospects for different target nuclei. Specifically, we compare 
three experiments with the same exposure time and target mass, but with xenon, germanium 
and argon nuclei (from left to right), in the planes defined by the lighter WIMP mass 
and scattering cross section (upper panels) and A versus again the scattering cross section 
(bottom panels). Although with naive assumptions on the background rate the heavier 
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the nucleus the better the discrimination prospects, we do not find any evidence that the 
discovery of a multi-partite dark matter is favored with one target nucleus choice versus 
another. The preferred mass and relative mass splitting for discrimination is also target- 
independent. 

IV. COLLIDER SEARCHES 

The new energy frontier that the Large Hadron Collider will soon explore will surlely 
have a tremendous impact on our understanding of physics beyond the Standard Model. As 
history teaches us, hadronic machines often correspond to "discovery" apparata, as opposed 
to leptonic colliders, where the goal has often been to attain precision measurements. In an 
optimistic scenario, physics results and discoveries with the LHC will trigger the interna- 
tional community to move on with a new generation of precision colliders, which many have 
envisioned as a natural next step, in the form of an international linear electron-positron 
collider (ILC). 

As discussed above, in the present model-independent setup, we specify neither the spe- 
cific particle physics model for the stable WIMPs nor the broader particle content of the 
model. The latter is of crucial importance to specify the decay channels and resulting 
observational signatures of new heavy particles produced at the LHC. Although studies 
have addressed this problem in the context of canonical collider signatures, including large 
transverse missing energy, jets and high energy leptons, where mass reconstruction of the 
produced particles can also be achieved (for a recent study see e.g. [27]), it is very problem- 
atic to set this problem in model-independent terms. We therefore focus on a simple scenario 
that an ILC might face: one where a massive charged particle species C is pair-produced 
in e + e~ collisions, and it decays into the two stable WIMPs with a certain branching ratio, 
giving off a high transverse-momentum lepton (for simplicity we assume the latter to be an 
electron or a positron). We showed schematic relevant Feynman diagrams in Figure [I] As 
for indirect detection, we assume there is no cross-channel where for instance the charged 
particle decays into both stable WIMPs, plus other standard model particles. This channel 
would anyway be kinematically suppressed. 

In this context, the key parameter of the game is the mass of the charged species 
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which will then determine the production cross section at the ILC 

r^oc-yi-— £. (20) 

Here we take the center of mass energy to be yfs = 500 GeV, as might be expected for the 
early phase of an ILC, and set a = a(mz) = 1/129. In the context of e.g. supersymmetry, 
the constant of proportionality in the above expression would be | q~, where is the charge 
of C. To be conservative, we take here such a constant to be 1 in our analysis. The additional 
parameters are the decay widths into the two stable species, plus light leptons (the case of 
quarks would be similar, in case C is a strongly interacting species): 



2 

m 



r L2 .. ,,. (l - ^ ] . (2i) 



As in the case of direct detection (and as discussed in Sec. II B), we can relate the to 
r a 2 ~ w, assuming that the particle that dominantly mediates the pair annihilation cross 
section into Standard Model particles is C (for instance, in the case of supersymmetry, a 
squark mediating the pair annihilation of neutralinos into a quark-antiquark pair in a t- 
channel). With the further simplifying assumption that C only decays into Xi or X2, we 
have the branching ratios 

BR lfl = =^=- . (22) 

J- 1 + 1 2 

Assuming that the cross sections are inversely proportional to the the dark matter energy 
densities we have a\ 2 oc <7i j2 oc l/f2i )2 , and, using (01), we can write 



BR\ = ^ 5- (23) 

1 + \/W § _ 2 



C 



The number of signal events is 



where 



BR 2 = 1 2 ¥ . (24) 

1 -I- 1 C I 
Jw m~—mZ 
v c 2 



N li2 = BR li2 N, (25) 



N = feS*- °c+C- X ^ ( 26 ) 

and we have the following estimates for the ILC parameters [28J 

f c s = 0.1 efficiency factor , (27) 

C = 100 fbT 1 integrated luminosity. (28) 
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We declare detection when either Ni > 20 or N 2 > 20. 

As for the discrimination of two distinct stable WIMP species in the final state, we 
consider the accuracy of ILC mass measurements, given the number of events and the 
intrinsic detector resolution: 

^ ~ Max ( ( ^) , *) (29) 
From [291 EOJ |3T] we get the following estimates: 

= 0.0005 (30) 

a = 0.11 (31) 



/ Ar»i, 2 \ 

V m l,2 / 



max 



where a indicates the statistically-limited detector mass resolution, as estimated in the 
simulations presented in Ref. [29]. Specifically, the simulations in [2S] present the relative 
error in the determination of the smuon mass for SPS point #1 [32J at a 500 GeV center of 
mass linear collider, including the effects of initial state radiation, beamstrahlung and beam 
energy spread. We then declare discrimination when both xi an d X2 are detectable and 

nil + 5Am! < m 2 — 5Am,2 (32) 

i.e. the two masses are discernible at 5cr. 

Fig. [9] illustrates the reach of an ILC with the present setup, in the (mi, mc) parameter 
space, for fixed values of the mass splitting A and relative dark matter species abundance w. 
While the blue edges essentially trace the kinematically defined parameter space region where 
rrtc > mi, the red contours correspond to satisfying both the requirement that mc > m 2 
and that a sufficiently large event statistics is achieved. In particular in our scenaro, for 
small w (for a Universe dominated by the heavier species) the heavier WIMP interacts more 
weakly (since it annihilates less efficiently) and C has a consequently smaller probability 



to produce it in the decay mode. This is also evident from Eq. (24) above. Comparing 
the upper two panels this is very easily seen: if w = 0.1 the charged particle must be light 
enough as to have a large enough number of events leading to the heavier WIMP in the 
final state (left panel). When w = 1, instead, the reach extends almost all the way to the 
kinematic limit. The other panels illustrate the same point: A smaller value of w (or a 
very large one too, see the bottom right panel with w = 10) indicates that one of the two 
stable species is the preferred decay mode for C, and it is harder to have enough statistics 

25 




250 



200 



O 



e 150 



100 




10 

m, (GeV) 



10 

m, (GeV) 



250 



£ 200 



E 150 



250 



200 



150 



100 



250 




m, (GeV) 



m, (GeV) 



FIG. 9: Collider searches, vs mi is plotted at fixed values of A and w. Color coding: blue stands 
for no detection, green for detection without discrimination, red for detection with discrimination. 
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FIG. 10: Collider searches, vs A is plotted at fixed values of m\ and w. Color coding: green 
stands for detection without discrimination, red for detection with discrimination. 



to achieve discrimination. If both particles are kinematically open final states, the w ~ 1 
scenario is the most favorable one at an ILC. 



Fig. 10 studies, for fixed values of the lighter WIMP mass (40 GeV and 99 GeV in 



the upper and lower panels, respectively) and for w = 0.1, 1 and 10, the regions of the 
(A,mc) parameter space where discrimination of two stable WIMPs is possible. Discovery 
is possible (green region) throughout the entire parameter space shown. The right (and 
lower, for mi = 99 GeV) edges of the red regions where the discrimination of two stable 
WIMP species is possible are dictated by requiring that the heavier WIMP be lighter than 
the charged particle, while those to the left depend upon the finite detector mass resolution: 
the best possible case (w ~ 1) suggests that an ILC might be able to discover two distinct 
stable WIMPs as long as their relative mass splitting is larger than 5-10%. The upper edges 
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COLLIDERS 




FIG. 11: Collider searches, vs w is plotted at fixed values of mi and A. Color coding: green 
stands for detection without discrimination, red for detection with discrimination. 

again are determined by the requirement of having a large enough statistics. The take-away 
message from Fig. [lOjis that the ideal relative mass splitting to single out two stable WIMPs 
at an ILC ranges between 0.1 < A < 2. 

We further study the role of w in Fig. [TTJ where we study discrimination of a dual- 
component stable WIMP scenario in the (w, mc) plane. The figure highlights that the ideal 
situation for discrimination corresponds to comparable decay branching ratios into the two 
stable species. The shape of the panels to the right, with larger A, is determined by the 
requirement of producing both WIMPs (e.g. in the lower middle panel, the second WIMP 
has a mass of around 200 GeV, hence the red region starts where C can actually decay into 
it). Again, when mc — > 250 GeV, the kinematic reach of the collider is reached, and the 
total number of events decreases accordingly, leading to worse discrimination performance. 
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FIG. 12: Collider searches. Here we replace A with m,2 in the parameter space (11) and plot rri2 



vs mi at fixed values of m^j and w. The asymmetry in the 4th plot comes from the combination 



of w and the masses in Eq. (23) and (24). Color coding: green stands for detection without 



discrimination, red for detection with discrimination. 



In the last figure of this section, Fig. [12j we study the plane denned by the two stable 
WIMP masses mi and m 2 , for fixed values of w and of mp. As expected, discrimination is 
problematic in the region of parameter space where mi ~ m 2 . In addition, the right edge of 
the plots to the right illustrates the critical issue of a largely asymmetric distribution for the 
relative dark matter densities. When dark matter is dominated by the first WIMP species 
(w = 10), we see from Eq. (23) that xi is only very rarely produced in the decay chain of 



C . The branching ratio BRi is even more suppressed if mi — > mg. Therefore, we see in the 
plots that at constant values of m 2 there is some upper limit to mi, beyond this limit we 
would not detect Xi anymore because the branching ratio BRi would be too small. 
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V. INDIRECT DETECTION 



The indirect detection of dark matter can proceed in a number of channels, or "mes- 
sengers" [33], including but not limited to antimatter, gamma rays and high-energy neutri- 
nos (for reviews in the context of several particle physics models for dark matter see e.g. 
[U [21 El Sj ) • The detection of gamma rays produced in the annihilation of dark matter, or in 
the radiation emitted in the subsequent energy losses of the energetic electrons and positrons 
(e.g. via Inverse Compton up-scattering of photons in the Cosmic Microwave Background) 
has been discussed at length in the literature (see e.g. the review of the possible dark matter 
search strategies with Fermi-LAT given in Ref. [31]). 

In order to discover a mult i- component dark matter scenario, we focus our attention on 
the detection of multiple monochromatic lines arising from the direct annihilation of dark 
matter pairs into two photons. This poses a fundamental problem. Additional lines can 
emerge from other channels, such as dark matter pair-annihilation into a photon plus a 
neutral gauge or Higgs Boson [35J, or in the context of universal extra-dimensional models 
[H [36]. Thus, given two lines for instance, one cannot distinguish whether the second 
line is due to a second WIMP or is just an additional line arising from the first WIMP 
(via annihilation into or H'y for example). Because of this caveat, indirect detection 
cannot be used to claim a discovery of mult i- component dark matter, but it could provide 
a valuable verification if evidence for such a scenario is found via direct detection and/or 
collider searches. Therefore, it is useful to carry out the analysis here and compare to the 
results obtained from the other channels in the previous sections. 

The signal flux of monochromatic photons for the XiXi ~^ 77 an d for the X2X2 77 can 
be cast as [37] 

_ 7.5 x 10~ 10 / w V ( (a XlXl _^ r/ v) \ /50 GeV\ 2 / AO" 



and 



cm 2 s \l + w J \10~ W cm 3 s -1 J \ m,\ ) \ sr 
7.5 xHT 10 w ( (a XlXl ^v) \ ( 50 GeV \ 2 { AVL 



X;i cm 2 s (1 + u;) 2 V10" 26 cm 3 s"V V( 1 + A ) m iy V sr 



J (AO) (33) 



J (AQ) , (34) 



where 



3 < ASi > S ( o.3 GeVen-O L ^ L ^ (Al)) ^ 

is the integrated line-of-sight dark matter density squared, over an angular region AQ, in 

a given direction in the sky. Given that the monochromatic two-photon line is typically 
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a statistically-limited indirect detection channel, the best channels include those involving 
observations with a large number of collected photons. Possible choices are the extragalactic 
diffuse signal from all redshifts and all halos, the galactic diffuse signal, the galactic center 
region and extragalactic targets such as local dwarf galaxies and nearby clusters of galaxies. 
We consider here study the Milky Way Galactic Center, which, for density profiles 

concentrated in the central regions of the Galaxy, can lead to the best signal-to-noise among 
the possible dark matter search targets listed above (see e.g. [38J). 

We assume here a fiducial NFW profile [39], featuring a local dark matter density at 
the Sun position of 0.3 GeVcm -3 , a scaling radius of r s = 20 kpc, and the high-energy 
Fermi-LAT angular resolution AQ = 3.83 x f0~ 3 sr, which corresponds to an angle of ip — 2 
degrees within the direction of the Galactic Center. We obtain the following: 

(~si~) J ~ 2,30 ^ AQ = 3,83 x 10 3 ^ 

As far as the background is concerned, we refer the reader to the detailed analysis of Ref . [38J ; 
here we take a subset of the background sources mentioned in that study, and consider contri- 
butions from (1) the H.E.S.S. detected source at the Galactic Center, (2) the EGRET source 
near the GC, (3) diffuse radiation, divided into a galactic and an extragalactic component. 
Specifically, we use the following background fluxes |i0] : 

/ f \ - 2 - 25 

fits = ID" 8 (^Vj GeV- cm- s~\ (37) 



egret = 2 - 2 X 10_7 (tt^) x ex P (-# 7 /30 GeV) GeV^ 1 cnT 2 s" 1 



VGeVy 

and, from [H] 



(38) 



<& kS ga i = H86 x 10- 6 (J^j GeV^ 1 cm" 2 s"\ (39) 
& - 3.66 x 10- (^f^) 21 ) GeV- cm" 2 8 - (40) 

We consider an effective area for Fermi-LAT of 8500 cm 2 , which is appropriate for the 
range of energy we are dealing with (see e.g. [12]), and an exposure time corresponding 
to 5 years with the GC within the Fermi-LAT field of view 50% of the time [13]. The 
product A e s x t Q bs is therefore « 6.7x 10 11 cm 2 s. Finally, we consider an energy resolution 
AE 1 /E 1 ~ 0.1. 
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Given a monochromatic flux of photons at a certain energy E , we compute the number 
of background events iVbckg in the bin between 0.95E and 1.05E , and compare it with the 



number of signal events iV S ig, computed from Eqs. (33) and (34). We then declare that a 
given dark matter setup is detectable if, for either xi or % 2 , the following two conditions are 
satisfied: 

Asi " > 5 (41) 



a/ -^bckg 

N sig > 50 (42) 

The two conditions above require a signal-to-noise (the latter assumed to be gaussian) at 
least at the 5 — a level, and a total number of signal photon events larger than 50. We 
instead declare to have discrimination for a setup where (1) both xi an d X2 are detectable 
and (2) the mass determination is such that 

mi + 5Ami < m 2 — 5Am 2 , (43) 

i.e. the two masses are descernible at 5cr. We estimate 

Am,- = — rrij (44) 
V 7 ^ 



Fig. 13 studies the (raj, (°~xrxi->~n v )) pl ane f° r fixed values of the mass splitting A and 
of the relative abundance of the two WIMPs, w. The dashed lines indicate the current con- 
straints on monochromatic gamma-ray lines as obtained in the analysis of Ref. [22] (see also 
[23]). and the regions above those lines are excluded. We are interested in looking for setups 
where the current exclusion limits fall well above the line that distinguish discrimination 
from detection, i.e. the boundary between the red and the green regions. 

First, we notice that in general, and not unexpectedly, the most favorable scenarios are 
those where w ~ 1, i.e. where the dark matter is roughly half made up of the first WIMP 
species \i an d half of the second species X2 (e.g. top three panels to the right). When the 
dark matter is mostly made up of one WIMP species, be it X2 (small values of w ~ 0.1) or 
Xi (large values of w ~ 10), discrimination is almost never possible. The one exception we 
find is the panel in the lower right corner, with A = 1 and w = 10, where discrimination 
would be possible with a pair annihilation rate within an order of magnitude of the current 
limits. Discrimination is also marginally possible for w = 0.1 and again A = 1 (second 
panel from the top in the left column). We thus conclude that for indirect detection, the 
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FIG. 13: Indirect detection. (a XlXl ^^v) vs mi is plotted at fixed values of A and w, as indicated 

in each panel. The current limits on gamma-ray lines are taken from the analysis of Ref. [22 . The 

arrows indicate the excluded region. Color coding: blue stands for no detection, green for detection 
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without discrimination, red for detection with discrimination. 



ideal setup for the relative mass splitting and abundance is when they both are close to 1, 
i.e. when the two WIMP species are relatively close in mass, and they contribute roughly 
equally to the dark matter content of the Universe. 
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FIG. 14: Indirect detection. {a XlXl ^^v) vs A is plotted at fixed values of m\ and w. The 
current limits shown are taken from [22J. The arrows indicate the excluded region. Color coding: 
blue stands for no detection, green for detection without discrimination, red for detection with 
discrimination. 



In Fig. [14] we see again that the best chances for discrimination of a dual-component 
dark matter scenario correspond to w = 1. In that case, the range of values of A for which 
discrimination is possible is rather wide, 0.01 < A < 10, but is optimal for A between 0.1 
and 1. As we found before, discrimination is also possible for w = 10 if A ~ 1 (see the 
bottom panel to the right in particular). 

In Fig. 15 we relax the assumption that (cr X2X2 ^ 77 w) = w x (a XlXl ^ yi v ), and study the 
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FIG. 15: Indirect detection. Here we relax the assumption (14) and we plot (o~ XlXl 



+ 77 



V) vs 



( a X2X2^77 v ) a ^ fixed values of m\ and A for equal energy densities U X1 = Q X2 or, equivalently, for 
w = 1. The current limits shown are taken from [22]. The arrows indicate the excluded region. 
Color coding: blue stands for no detection, green for detection without discrimination, red for 
detection with discrimination. 



((cr XlXl _^ 77 t> ), x (c X2X2 _^. 77 u)) plane. We fix m\ = 40 GeV in the upper panels and 100 GeV 
in the lower panels, and set A = 0.1, 1 and 10 from left to right. In all cases, we find 
that there is a range of values for the cross sections where discrimination is possible and 
consistent with current constraints. This range typically involves cross sections of at least a 
fewxlO -26 cm 3 /s ballpark, and extends up to rather large values (particularly when one of 
the two species is subdominant in the overall dark matter density). We therefore find that 
for WIMPs with total pair annihilation cross sections close to the typical values needed to 
achieve thermal production in the early Universe (i.e. (erf) ~ 3 x 10~ 26 ) the branching ratio 
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into the 77 channel must be extremely large, if not the dominant one. 



In Fig. 16, we study the plane defined by the two WIMP masses (we relax here the 
assumption that m 2 > mi), and assume a fixed pair annihilation cross section (cr XlXl _, 77 t>) = 
1CT 26 cm 3 /s. We then simply consider four different values of w, a quantity which also sets 
the pair annihilation rate into two photons of the second WIMP. While in the top left 
panel the dark matter is almost entirely made up of X2, and that's the only species that 
can be detected in the channel we consider here (in the bottom right the same is true for 
Xi), the most promising cases are again the intermediate ones. For the particular choice of 
parameters we made, w = 10 does not lead to any viable discrimination parameter space 
which is not already excluded by current data, while the "symmetric" w — 1 case allows 
for discrimination as long as both WIMPs are (1) not degenerate (green line at mi = m 2 , 
where clearly no discrimination is possible) and (2) both lighter than roughly 15 GeV. This 
precise mass range obviously depends upon the assumed value for (a XlXl ^ 77 t> ). 

This point is shown explicitly in the bottom four panels, where we set (cr XlX1 _, 77 u) = 
10~ 25 cm 3 /s (top) and (o" XlXl _» 77 f ) = 1CT 24 cm 3 /s (bottom). Again, we find that the ideal 
situation is a symmetric case (w = 1), with masses that can be detectable up to 150 GeV for 
a cross section of 10™ 25 cm 3 /s and above a TeV for 10~ 24 cm 3 /s. With w = 10 discrimination 
is possible only for large masses and for very large pair annihilation cross sections into two 
photons (10~ 24 cm 3 /s, bottom right panel). 



VI. RELATED WORK AND POTENTIAL MODELS 

The possibility of a "complex", multi-partite dark matter sector has been envisioned 
in different contexts, including models where the dark matter is partly made of WIMPs 
and partly made of e.g. axions [44] or sterile neutrinos [15]. In addition, models with 
multiple WIMP species as dark matter candidates have been discussed, motivated either by 
theoretical considerations (see e.g. (HJl EH SHI EH]), or to explain anomalous signatures in 
cosmic ray or direct detection experiments that could not be easily accounted for with a 
single dark matter particle (e.g. [501 ED 152]). 

Without claiming to be exhaustive in our account, we summarize here a few recent studies 
that addressed, in specific particle dark matter models, the possibility that the dark matter 
sector is comprised of more than one, stable weakly interacting massive particle (we therefore 
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do not count here as "multi-partite" (WIMP) dark matter models those frameworks where 
e.g. besides a stable WIMP one has a contribution to the dark matter density from neutrinos 
(for instance in mixed warm-cold dark matter models with sterile neutrinos jl5]) or from 
axions [S]). 

Ref. [50J showed that a two-component dark matter scenario (invoked in their analysis 
to explain both the Integral/SPI 511 keV line in the galactic bulge and the EGRET GeV 
excess, the first anomaly requiring a light dark matter WIMP in the MeV range and the 
second one in the GeV range) can emerge in the context of N = 2 supersymmetry with two 
discrete symmetries, the "usual" R parity of the MSSM and an additional M-parity. The 
lightest R— and M— parity odd particles are stable, the latter in principle have desirable 
features to fit the Integral/SPI 511 keV line. 

In Ref. [IB], a scenario was proposed for the radiative generation of Majorana neutrino 
masses via the introduction of additional superfields to the MSSM, namely three heavy 
neutral Majorana fermion singlets, two more SU(2) Higgs doublets (in addition to the usual 
two MSSM Higgs doublets), one scalar singlet. These superfields are odd under a new 
discrete parity Z' 2 . The model not only generates naturally small neutrino masses, but it 
can also account for the baryon asymmetry via leptogenesis, and provides a multi-partite 
dark matter scenario involving the lightest odd-charged particles under the two discrete 
symmetries. The phenomenology of this scenario, as well as of the more general setup 
where the dark matter is comprised of a fermion singlet and of a scalar singlet, with quartic 
interactions with the SU(2) Standard Model Higgs, was studied in detail in Ref. [H]. That 
study makes several specific and model-dependent assumptions on the particle physics model, 
and uses a different parameterization and definition of whether a model is or not detectable. 
A direct comparison of the results presented here with those in Ref. [IB] is thus arduous, 
although Ref. [IE] presents a specific model realization of the setup we hereby investigated. 

A scenario that incorporates a solution to the problem of small neutrino masses, low- 
scale leptogenesis and a multi-partite dark matter scenario, involving a Z2 x Z' 2 discrete 
symmetry, was also recently envisioned in Ref. [53]. The two dark matter particles of that 
scenario are a gauge-singlet Majorana neutrino S and a singlet scalar <ft. In addition to 
the phenomenology of the model at direct dark matter detection experiments and with 
colliders (specifically, through the search for invisible Higgs decays at the LHC), Ref. [53] 
also discusses the possibility - which we do not entertain in the present analysis - of pair- 
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annihilation processes of the type SS — > 4>4> or , depending on the mass spectrum, 00 — > SS. 
Those processes can play an important role in setting the relic abundance of the two dark 
matter particle species in the early universe [53] . 

An interesting realization of a two-component dark matter model was presented in 
Ref. [3?] , which focused on an extension of the MSSM featuring a spontaneously broken 
U(l)' gauge symmetry giving rise to a conserved residual discrete [/-parity, and to an addi- 
tional stable WIMP that can contribute to the dark matter in addition to the ordinary LSP. 
This possibility is soundly motivated, e.g. from possible solutions to the MSSM \x problem 
and to the suppression of higher dimensional i?-parity conserving operators that can medi- 
ate lepton and baryon number violating interactions. Ref. jl?] considers the phenomenology 
of direct detection and collider searches for both the case where the scalar or the fermion 
component of the [/-parity odd superfield is the additional lightest stable particle. 

Ref. [49J pointed out that in gauge mediated supersymmetry breaking it is natural to 
expect multiple stable hidden sector particles that quite remarkably have the desired thermal 
relic abundance to explain the dark matter density in the Universe. Along similar lines, this 
possibility was further investigate in Ref. [53]. The recent lepton anomalies reported by 
the Pamela [55J and ATIC experiments (the latter only seen to a lesser extent by the 
Fermi gamma-ray space telescope, (S3 ES]) have triggered model-building that included the 
possibility of a multi-partite dark matter setup. This includes e.g. the study of Ref. [5T] . 
where the large boost factors needed for the annihilation rate of dark matter particles needed 
to explain the lepton anomalies are made compatible with a large enough relic abundance 
via the decays of a meta-stable heavier particle species, that is the dominant dark matter 
component prior to its decay. Ref. |52j studied instead the possibility that several anomalous 
signals that might be explained with two particle dark matter particles: one is a heavy 
quasi-sterile neutrino, pair annihilating preferentially into leptons, as needed to explain the 
mentioned lepton anomalies, and the second one is a possibly lighter particle that could 
account for the DAMA signal [511 [521159]. 

Unlike in the work summarized above, we have taken in this study an entirely model- 
independent approach, and do not refer to any specific multi-partite particle dark matter 
setup. To a good approximation our findings can be appropriately applied and translated 
to any of these particular models. 
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VII. CONCLUSIONS 



We have focused in this study on the prospects for establishing, with dark matter searches 
and a future linear collider, the existence of dark-matter physics featuring two stable weakly 
interacting massive particles in our Universe. Here we summarize our general results, and 
outline, under our assumptions, the generic requirements on the parameter space of this 
scenario in order for significant experimental signatures to be imprinted upon one of the 
channels under consideration. 

• Direct Detection: we considered the possibility of discerning a multi-component dark 
matter scenario via fits to recoil spectra in direct detection experiments with a setup 
of 5 yr of data taking with a 1000 kg target. We found that: 

— To have discrimination of two dark matter species at least one of the WIMPs must 
have a light mass, in the 10 GeV to 30 GeV range, even with large detectors. 
The mass of the heavier species can be as large as a TeV, for sufficiently large 
experimental setups and scattering cross sections. 

— A necessary condition for the discovery of multi-partite dark matter with direct 
detection experiments is that the two species contribute comparably to the overall 
dark matter density of the Universe. 

— The relative mass splitting A for which we find that it would be possible to dis- 
cover multi-component dark matter is rather narrow for direct detection, ranging 
between 0.5 and 2, and peaking where the mass ratio between the two dark matter 
species is in the proportion of 2:1 (A = 1). 

— We do not find the choice of the nuclear target to be crucial to discover a multi- 
component dark matter scenario. 

• Linear Collider, we assumed that in a future electron-positron linear collider with a 
center of mass energy of 0.5 TeV a pair of charged massive particles is produced and 
decays into one of the two stable WIMPs plus a light lepton, with a branching fraction 
related to the relative WIMP abundance in the overall dark matter density. We found 
that: 

— When the decay of the massive charged species into the WIMPs is kinematically 
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open for both, the discovery of two stable species (on collider scales) is possible 
for relative mass splitting A as small as 0.05 and up to 2. 

— The best case scenario is again when the relative abundance of the two species is 
comparable, which, with our assumptions, corresponds to comparable branching 
ratio into the two WIMPs for the decay modes of the charged massive species. 

— Discovery prospects worsen near the kinematic limits of the collider, simply due 
to lack of statistics, but the reach extends almost all the way to the kinematic 
limit in most cases. 

• Indirect Detection: we considered the detectability of two distinct gamma-ray lines 
associated with the direct annihilation of the WIMPs into monochromatic gamma 
rays. Because of the inability to distinguish whether the second line comes from a 
second WIMP or from a different decay mode of the first WIMP, this search channel 
cannot lead to a discovery in its own right. Nevertheless, it could prove useful for a 
verification of the dual-component dark matter scenario. We found that: 

— Large pair annihilation cross sections in the 77 channel are needed to have a de- 
tectable signature of two stable WIMPs, typically of the order of 10~ 26 cm 3 / s, or 
at most two orders of magnitude below the current, pre-Fermi-LAT experimental 
sensitivity. 

— Under the assumption (which could be relaxed) that the total annihilation cross 
section is proportional to the annihilation cross section into 77, this detection 
channel strongly depends upon the relative composition of the dark matter; the 
possibility to find two stable WIMP species is strongly favored if the two species 
contribute almost equally to the overall dark matter density. 

— Indirect detection of dual-component dark matter is quite effective in a wide 
range of masses, and for a considerably large range of relative mass splittings 
between the two species' masses, A, which can range between 0.01 up to more 
than 10. 
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